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Abstract 



We describe a formalism for solving the BFKL equation with a coupling that runs for 
momenta above a certain infrared cutoff. By suitably choosing matching conditions 
proper account is taken of the fact that the BFKL diffusion implies that the solution 
in the infrared (fixed coupling) regime depends upon the solution in the ultraviolet 
(running coupling) regime and vice versa. Expanding the BFKL kernel to a given 
order in the ratio of the transverse momenta allows arbitrary accuracy to be achieved. 
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Diffractive processes are described, at least within the framework of perturbative 
QCD, by the BFKL equation This equation sums the leading and next-to-leading 
Ins [fj] terms for the amplitude involving the exchange in the t— channel with the 
quantum numbers of the vacuum. 

In the case of fixed coupling, the energy dependence of these diffractive processes is 
determined by the emergence of a branch cut in the Mellin transform plane (u— plane), 
with a branch point at u > 0. Lipatov || has demonstrated that accounting for the 
running of the strong coupling together with some (non-perturbative) information 
about the infrared behaviour of QCD leads to a separation of the a;— plane singularity 
structure into a series of isolated poles. More recently, Thorne H] has shown that if 
one uses the solution to the BFKL equation with running coupling, the factor in the 
amplitude which is totally calculable in perturbative QCD has an essential singularity 
at u — 0, such that there is no power-like energy dependence. This is in accord with 
the recent work of Ciafaloni, Colferai and Salam, who have shown that the leading 
cu-plane singularity is of non-perturbative origin ||. 

The BFKL equation is a diffusion equation and so there will always be some 
contribution from the infrared region. One way of treating this region is to assume 
that the strong coupling ceases to run below a certain scale. In this case, we expect 
that for sufficiently low values of the transverse momenta of the exchanged gluons, 
the s— behaviour will be dominated by the fixed coupling solution, whereas for large 
enough transverse momenta the solution of Q will dominate. We should like to 
understand the smooth extrapolation between these two extremes. 

In this letter we describe a formalism for solving the BFKL equation with strong 
coupling running above some infrared scale, in which we can see how the transition 
from the fixed coupling to the running coupling solution arises. 

Formalism for Solving the BFKL Equation with Running Coupling 

To begin with we consider only the leading order BFKL kernel and describe later 
how the formalism may be extended in order to include the substantial part of the 
higher order corrections. 

Since the BFKL equation is a diffusion equation there is always diffusion into the 
infrared regime, where renormalization group improved perturbation theory breaks 
down and one needs supplementary information about the behaviour of QCD beyond 
perturbation theory. For our purposes we simulate this infrared behaviour by assum- 
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ing that the coupling freezes below a critical value of t = \n(k 2 / Aq CD ), t = t . We 
also make the assumption that a s (t ) is sufficiently small that perturbative results in 
the infrared regime maintain some level of credibility. 

The full BFKL equation is solved in terms of a complete set of eigenfunctions 
fu){t) of the kernel with eigenvalue u. For any oo there will be some diffusion into the 
infrared regime. We therefore write f w {t) as 



and the kernel as 
K(t,t') = 



U(t) = f<(t)9(t -t) + f>(t)6(t-t ) 



^/C^MWo-O + ^WW'-^o)] 0(t-t') + 



(1) 



^-/C>(t, 9(t - t') + ^ 7 /C>(t, f) 6(t> - t ) 



9(t' - 1) 



(2) 



where b = Po/12, (3q = 11 — 2n//3, being the first coefficient of the f3— function for 
the running of the QCD coupling. We define /C < (t,t') and /C > (t, t') as the limits of 
sums, i.e. 

/C<(M') = lim /C<(M') 



N->oc 

K>(t,l!) = lim /C^(M') 



N-+oc 



where 



^(t,f) = El^^-^it-t')} 

r=0 L r + 1 J 

N ( 1 1 

^(MO = eK +1/2)( ^' ) --xt^-o}. 



(3) 



(4) 



The functions e 7 * are eigenfunctions of these kernels with eigenvalues Xn(i) an d 
Xn(i) respectively, i.e. 

ft 



/C<(M')e 7i V = xUl^ 



where 



N 



TV ( x 



^[r + 1/2 + 7 r + 1 

M 

2-7 r + lj 



(5) 



(6) 



(7) 



Note that X©i(7) = 7e — VKV^ + 7) an d XSo(7) — 7b — ^(1/2 — 7), and we have the 
leading order BFKL kernel. 
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In this paper we truncate the kernel after N terms. Note that the case N = 
reduces to the collinear model of ||. For this truncated kernel the spectrum of 
eigenvalues is 

N n N 2 

-Y^— <bt co <Y- -. (8) 

+ 1 ~ _ ^(2r + l)(r + l) 1 ' 

The infrared part of the eigenfunction f^(t) obeys the integral equation 
bt uf<(t) = f K,<{t,t')f<{t')dt' + I*" K>{t,t')f<{t')dt' 

J — oo J t 

+ r mt^fiwdt' (9) 

whereas the ultraviolet part of the solution f^(t), obeys the integral equation 

ft roo f 

btuf>(t) = / K<(t,t')f>(t')dt' + / K>(t,t')-f>(t')dt' 

Jt Jt t 

+ r ic<(t,t')f<(t')dt'. (io) 



We note that these equations are both inhomogeneous, reflecting the fact that even 
for t < to, where the coupling is fixed, the BFKL equation involves diffusion into the 
running coupling regime and vice versa. 

We can convert these two equations into homogeneous integro-differential equa- 
tions by operating on both sides of eq.([|) with the operator 

and operating on both sides of eq. ( |I0|) with the operator 

°?< f > = l(H-(5 + JV ) e ' /2 

These equations are (N+l) th order in derivatives and so each one of these has solutions 
which contain iV + 1 arbitrary constants of integration. However, these constants are 
fixed by the requirement that the original integral equations (0,0) mus t be obeyed, 
together with the requirements the eigenfunctions should be normalizable, which 
means that f^(t) must be square integrable as t — ► — oo and f^(t) must be square 
integrable as t — > oo. 

The general solution to the integro-differential equation for f^(t) is 

N+2 
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where 7° are the solutions to 

but = Xn(i) + Xn(i) = Xn(i) (12) 

with 3?e(7) > to ensure square integrability as t — > — 00. In the allowed range of 
eigenvalues, eq.(|8|), there will be iV + 2 such solutions. For example, in the case N = 
the solutions are 

7? = ±*\/(V£) (13) 

where 

A„ 1 



but + 2 

For the allowed range, A > 1/4, and so these values are purely imaginary. In the 
case N = 1 we have 



7? = ±\/ (y " 2Ai ± y/l - 2\ x + 4A? ) (14) 

where 

Ai 



but + 3 

For the allowed range, Ai > 3/16 and we see that two of the solutions are purely 
imaginary, Hfx, and two purely real ±/?i. Thus the general solution is 

/<(*) = axe^ + a 2 e^ + a 3 e^. (15) 

For the integro-differential equation for f^(t), it is convenient to work in terms of 
the Laplace transform ^(7) defined by 

nr " as/*'-*** (16) 

in which the integral over 7 is performed over some suitable contour, described below. 
The function f w (j) obeys the differential equation 

w6/J(7) + XN^fM+X^Ul) = 0, (17) 

where the prime indicates differentiation with respect to 7. This equation can be 
solved in WKB approximation yielding the result 

L(i) = ^^ e ^-i/V y)+ ^^M- (l8) 

where 

V( 7 ,u) = xUi)-^Xn'{i)- 
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Note that if we set the term Xn'(i) to zero we recover the expression obtained in M 
for the case of deep-inelastic scattering in which the transverse momenta are ordered 
and the coupling is chosen to run always with t. 

The solution fllTf ) is a valid approximation except near the points where V{^) 
vanishes. We could improve the solution by expanding V{^f) to linear order around 
the zeroes and matching the solutions either side to an Airy function. However, the 
singularities in f w at these zeroes are integrable so that such an improvement would 
have negligible effect on the inverse Laplace transform. 

For t < t c the solution is oscillatory. The value of this critical value t c can be 
estimated by approximating the integral of eq. (|l6l) by the saddle point at 

2ubt = X n(i) + y/vHu). 

t c is the minimum value of t for which this has a solution for purely real 7. Again, if 
we were to neglect Xn'(i)-> this would occur at 7 = and we have approximately 

buo 

This oscillatory part is "matched" to the fixed-coupling solution by requiring 
consistency with the original integral equations flj^, [10]). To see how this is possible, 
we need to examine in more detail the possible contours of integration for eq. flT6|) . 

The 7-contour 

From the positions of the simple poles in Xn{i) an d Xn(i)i we n °t e that /a, (7) 
has branch points on the real axis at 7 = r + 1/2 and 7 = — 1/2 — r for r = 0, ■ ■ ■ N. 
The branch cuts can be "combed" in the direction of the positive or negative real 
axis in such a way as to leave a portion of the real axis between n — 1/2 and n+ 1/2 
for which foil) is analytic, n runs from —N to N. 

For each such "combing" a contour, C n can be chosen that crosses the real axis 
between n — 1/2 and n + 1/2. Since it > we require that the ends of the contour 
turn over so that at the ends of the contour 9fte7 — >• —00, thus ensuring that the 
inverse Laplace transform integral exists. The integral over each such contour is a 
valid solution to the integro-differential equation for f>(t) and since this equation is 
linear, the most general solution which is square integrable as t — > 00 is the sum of 
the integrals over these contours with —N < n < 0, with arbitrary coefficients b n , 
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i.e. 



n=-N 



E b n / ^7^(7)e 7t , (19) 



with fu){l) given by eq. (|T8|) .f| As an example we consider the N = 1 case. The four 
"combings" are shown in Fig. p] together with the four contours, C n . Only C_i and 
C give a square integrable solution, so the general solution is 

{/■oo _, /"OO _ ~| 

6-1 / dve- 1 j^iv-iy^ + 6 / dvUiv)e M . (20) 
J -co J —oo ) 

Returning now to the original equations (§) and (0), we see that these are obeyed 
by the solutions ( |19D and ([TT|) provided 



N+2 7 ?to f C^ e 7*o 

S> (7 ._ r _ 1/2) -'° X>/ c . - ° <'=<>•••">■ < 21 > 



iV+2 7 > /"'Me 7 * 

^ 3 (l° + r + l/2) n t^ N n Jc n y ( 7 + r + l/2) 1 J 1 ' 

Thus we obtain 2N + 2 relations between the N + 2 coefficients cij and the N + 1 
coefficients (the overall normalization of these coefficients is determined by the 
normalization condition of the eigenfunctions). This then determines uniquely the 
required eigenfunctions. 



Equations (g, [L0|) guarantee the continuity of the eigenfunctions at t = t , i.e. 
/of (^o) = /ilT (^0)5 but the derivative is not continuous. This results from the fact that 
we have taken an expression for the coupling as a function of t, whose derivative is 
not continuous at t = to- 

We can now see how the inverse Mellin transform, f(s,t) of the solution f w (t) 
extrapolates smoothly between the power-like behaviour for t ~ to and the softer 
behaviour (with at most logarithmic dependence on s) for t ^> to- 
Examining the 7— integral for one particular contour, C n , and for t > to, we see 
that this will be dominated by a saddle point j s (u,t), which is the solution to 



2b ut = XN^ + y/Vi^u). (23) 

This is multi-valued and the dominant saddle point will be the one closest to the 

integration contour. We call this saddle point Js(uj, t). In other words, f^(t) of 

tWe have absorbed the factor of 2-rri in the denominator in eq.([l^) into the definition of the 
coefficients b n . 
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Figure 1: The four contours for the case N = 1 corresponding to the four different 
ways of "combing" the branch cuts. 



eq.flTED from the contour C n will contain an oj— dependent prefactor, which may be 
written (after a suitable integration by parts) 



The essential singularity at oj = is now encoded in j™(u,t) which acquires a singu- 
larity as oj — > 0, i.e. the R.H.S. of (p3|) can only tend to zero as |7| — > 00. 

Now examining eqs.(^l|, |2^) , we see that the contour integrals multiplying the 
coefficients b n are also dominated by a similar saddle point but with t set to to- Thus 
these coefficients appear in eqs.(|H|, with a prefactor 

exp( [ t0 ^(aj,t')dt' 



Therefore these coefficients also possess an essential singularity at oj = 0. Furthermore 
b n will contain factors of the form exp(7o(cu)to) from the L.H.S. of ( pi , f22|) . 



For t close to to each term in the sum of eq.(|T^), i.e. the product of b n and the 
contribution to f^(t) from the integral over the contour C n will have a dominant 
oj— dependence of the form 



exp 



{7 s >,t )(t-to)+76V)to}. (24) 



Note that due to the partial cancellation, the coefficient of the term that gives rise to 
the essential singularity vanishes as t — ► t . 

For sufficiently large s, the inverse Mellin transform probes small oj, where the 
term in the exponent proportional to 7" to) will dominate even if (t—to) is small and 
the s-behaviour is dominated by the essential singularity in the Laplace transform, 
fuiij)- However, for more moderate values of s we can neglect this term when t is 
sufficiently close to to an d the s— dependence will be dominated by the fixed coupling 
power-like behaviour. As t is increased away from to, the boundary in s where the 
soft behaviour takes over from the power-like behaviour decreases and we see that as 
t — > 00 we recover the solution of J|. 

This formalism is considerably simplified if we assume that to is sufficiently large 
that exp(— 1 ) <C f . In this case we only need consider the right-most contour consis- 
tent with square-integrability as t — ► 00. This is the contour Cq. In this case, within 
the approximation of performing the 7— contour integral and the inverse Mellin trans- 
form by the saddle-point method, the s— dependence may be written as s w ( s '*), where 
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Figure 2: Q(s,t) against t for Ins = 23 

u(s,t) is the solution to 

ln(s) = -^^(t-to) j\° s (u,tyt' + ^(uj)Mm(t,t )Y (25) 

and 7™(cu) is the solution to eq. fll2]) with the largest real part. In Fig. || we plot 
u(s,t) against t for a large value of s (Ins = 23). For simplicity we have taken the 
N = case, although we expect the result to be qualitatively the same for higher 
values of N, and we take to — 4. We see that it is almost constant for t < to and 
diminishes as one increases t from to, indicating the smooth transition from hard to 
soft s— dependence as t is increased. The discontinuity of the slope at t = t is a 
reflection of the fact that we have taken a sharp discontinuity in the slope of the 
running coupling at to an d would be smoothed out if a smoother transition were 
taken. 

The asymptotic behaviour is only obtained for much larger values of t and in Fig. 
|3] we show the s dependence of u)(s,t) for t = 40 and notice that it is much closer 
to zero than the case where t ~ to and decreases as In s increases, resulting in an 
s— dependence for the amplitude that has no fixed power of s. 

Accounting for Higher Order Corrections 

The higher order BFKL kernel does not lend itself readily to an expansion in 
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Figure 3: u(s,t) against Ins for t = 40 

powers of e*~*'. Furthermore, it is known that these corrections are large. 

However, it has been pointed out || that a substantial part of the higher order 
corrections consist of "collinear corrections" which are required to ensure correct 
behaviour when the formalism is applied to deep inelastic scattering. It was shown in 
H and confirmed in |7J] that once these corrections are accounted for, the remaining 
higher order corrections are modest. The simplest way of encoding these collinear 
corrections and summing them is to replace 7 by (7 — oj/2) in Xn{i) an d by (7 + u;/2) 
in X >( 7 ) g|. 

This formalism is readily adapted to the procedure described here. In the expres- 
sion ( |Tl~D for fu(t) the quantities 7°(co>) are taken to be the solutions to the implicit 
equation 

X^( 7 -^/2) + x^(7 + ^/2) = cobto, (26) 
and the expression for the Laplace transform of f> (t) becomes 

~ Ul) = WWW 11 ^ f ( X ^ 7 ' ~ Uf2) + X>n(1 ' + ^ /2) + V / W^)) > 

(27) 

where 

V{i,u) = ( X ^(7-^/2) + Xv(7 + ^/2)) 2 -4^x^(7 + ^/2). (28) 
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Provided that we restrict our analysis to the case < 1, the singularity structure 
and consequently the analysis of the contours of integration remains unchanged from 
the leading order case, although the exact positions of the branch points will have 
moved. 
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